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Exclusive electroproduction of lepton pairs as a probe of nucleon structure 
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We suggest the measurement of exclusive electroproduction of lepton pairs as a tool to study 
inter-parton correlations in the nucleon via generalized parton distributions in the kinematical region 
where this process is light-cone dominated. We demonstrate how the single beam-spin asymmetry 
allows to perform such kind of analysis and give a number of predictions for several experimental 
setups. We comment on other observables which allow for a clean separation of different species of 
generalized parton distributions. 
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Since the discovery of a non-zero spatial extent of the 
proton in pioneering measurements of electromagnetic 
form factors pL the exploration of the hadron's inter- 
nal structure in terms of quark-gluon degrees of free- 
dom was the subject of intensive theoretical and exper- 
imental studies. Until recently, however, only informa- 
tion on one-dimensional slices of the nucleon was ex- 
tractable from hadronic observables: The spatial charge 
and magnetization distributions from the form factors 
alluded to the above; the a;-momentum fraction space 
parton densities from inclusive reactions, e.g., deeply in- 
elastic scattering of leptons off hadrons g. With the 
realization of the power of exclusive reactions B, the 
opportunity of nucleon holography [Q was brought to 
life; namely, the simultaneous measurement of the lon- 
gitudinal momentum fraction of partons and their trans- 
verse localization within the area of resolution set by 
the photon virtuality || (see also ||, 0])- The corre- 
sponding hadronic characteristics, known as generalized 
parton distributions (GPDs), intertwine the aforemen- 
tioned conventional hadronic functions. However, they 
also depend on an extra variable: the ^-channel longi- 
tudinal momentum fraction r), the so-called skewedness, 
which changes the apportion of longitudinal momentum 
between the absorbed, (x — r]), and created, (x + rj), 
partons. Thus, GPDs, contrary to conventional Feyn- 
man's parton densities, give access to inter-parton corre- 
lations rather than mere probabilities through the study 
of skewedness dependence. 

The pioneering experimental study of these new func- 
tions at HERMES §, CLAS § and HERA |o| has 
been done through the exclusive leptoproduction of the 
real photon off the proton Ep — > E'jp' at high momen- 
tum transfer. This process is sensitive to deeply vir- 
tual Compton scattering (DVCS) amplitude (see Eq. (B)) 
parametrized in terms of GPDs. The latter enter also in 
the leptoproduction of mesons, i.e., Ip — > E'Mp'. Sin- 
gle lepton and nucleon spin asymmetries O, [l2| in the 
former case allow the direct measurement of GPDs, how- 
ever, in a very specific kinematics, namely, x — —rj. This 
restriction is a direct consequence of the reality of the 
outgoing photon. Such a measurement is unable to con- 



strain the spin sum rule involving the parton's orbital mo- 
mentum jn|. Here a processes is needed where the mo- 
mentum fraction varies independently of the skewedness. 
This can be achieved by the relaxation of the reality 
condition for the photon in the final state, i.e., by the 
measurement of lepton pairs produced from a timclike 
7-quantum in the elastic electron scattering ep — > e'p' EE 
p"2| , |i"4j ; see also Ref. ITJ] for a related reaction. 




FIG. 1: Fhe kinematics of the lepton pair production 
in elastic electron scattering off the proton, e(k)p(pi) — -> 
e{k')p(p2) £{£-)£{£+), in the fixed target setup. 

The differential cross section for the process in the 
kinematics displayed in Fig. [l] reads 



da = 



16(2tt) 3 Q 2 



lAfdxBdyd^ld^dM^dfii. (I) 



The incoming photon virtuality and the invariant mass 
of the lepton pair are q\ = — Q 2 and M 2 - t = (E + + EJ) 2 , 
respectively. The Bjorken variable and the lepton energy 
loss is defined conventionally as x B = Q 2 /(2pi ■ qi) and 
y = p ■ qi/p ■ k. A 2 = (jp2 — Pi) 2 is the i-channel mo- 
mentum transfer. The solid angle of the lepton pair in 
the photon rest frame is dflg = sinOfdOfdip^. (j) is the az- 
imuthal angle between the lepton and hadron scattering 
planes and plays a distinguished role below. The ampli- 
tude A represents the sum of the signal in question and 
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FIG. 2: Contributions to lepton pair production: (from left 
to right) virtual Compton and Bethe-Heitler processes. 

"contaminating" Bcthc-Heitler processes, 

A = Avcs + -Abhi + Abh 2 , (2) 

corresponding to the first, second and third diagrams in 
Fig. 0, respectively. The latter two are expressed in terms 



of the Dirac and Pauli form factors parametrizing the nu- 
cleon matrix element of the quark electromagnetic cur- 
rent, 

J„ = U2 (Vi + ^g^F?) (A 2 )-i , (3) 

with Ui being the nucleon bispinor itj = u(pi). As we 
already mentioned above, they are measurable elsewhere, 
see, e.g., [|l6| for the most recent results. 

The gauge invariant decomposition of the hadronic ten- 
sor contributing to -4vcs was found in Ref. jl7| by an 
explicit twist-three analysis at leading order of perturba- 
tion theory. To leading twist accuracy, it reads 
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Here the symmetrized combinations of momenta q = 
(qi + Q2 ) /2 and p = p\ +p 2 are used to define the gener- 
alized Bjorkcn variable £ and skewedness parameter r\: 
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Since the light-cone dominance in Eq. (j^j is set by the 
average virtuality q 2 at small A 2 and moderate x B , the 
perturbative QCD approach is applicable provided 



q 2 \ = \Q 2 - M'ff + A z /2| 12 > 1 GeV^ 



(7) 



It is lost for Q 2 ~ which implies, according to Eq. 
(|J), that £ approaches zero. The minimal allowed value 
is set by |£ m in| ~ 1 GeV 2 /ys, where ^/s is the center-of- 
mass energy. In the light-cone dominated region ^ , the 
Compton form factors factorize into calculable coefficient 
functions and GPDs. To leading order in coupling con- 
stant the Compton form factors satisfy the generalized 
Callan-Gross relation 

V 2 = £Vi . (8) 
Performing the Dirac decomposition one gets 

-E) (x, v ,A 2 ) Ul , 



Vi=u 2 dxC.(x,0[l+H 
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in terms of GPDs H, E, H and E ||. The plus subscript 
stands for the contraction of the corresponding Lorentz 
index with the light- like vector = —£(2^ + £p M )/<7 2 , 
which projects out the leading power contribution. From 
here, one immediately sees the difficulty to measure these 
functions: one of the dynamical variables enters inte- 
grated out with the coefficient function, which reads to 
leading order in QCD coupling constant, 



C T (x,0 



1 



£-x~iO £ + x-iO 



(9) 



One can get rid of the convolution provided the observ- 
able is sensitive to the imaginary part of the Compton 
form factors only, i.e., 3mVi, 3m.4i. 

The most illuminating experimental observables in this 
respect are single beam or target spin asymmetries. Their 
advantage is that they (i) depend linearly on Compton 
form factors, and (ii) are proportional to their imagi- 
nary part. The complete result for the cross section 
([!]), which is represented by a double Fourier sum in 
azimuthal angles da ~ J2 m n cos ( m 4>){ cc m,n cos(m<pf ) + 
cs m ,n sm(mipe )} + (cos «-> sin), will be published else- 
where. Here we limit ourselves to the case of the lep- 
ton helicity difference and when one does not distin- 
guish the final state pair with respect to their angular 
distribution. The integration over the solid angle 
leads to the vanishing of the interference term AvcsA Bii ^ 
and so the only term which survives in the asymmetry 
Act = <t(A = 1) - o{\ = -1 ) is ^vcs^b Hi : 
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dAa 2ai m ( 2 - v)V\l^ - x b)Q 2 - xbM^tJA^ - A 2 sin 



dxBdyd\A 2 \d^dM a 3tt Q 2 M 2_ A 2 (^/I" 17(Q 2 + + 2(2 - y) ^(1 - xb)Q 2 - x B M 2 ^A 2 min - A* cos 

x (^F 1 (A 2 )^(e,ry,A 2 )+^i+^)(A 2 )iJ s (e, ? ?,A 2 )-^^(A 2 )ii; s (e,r,,A 2 ^ , (10) 



where F S (£, V ,A 2 ) = F(£,r?,A 2 ) - F(-£,77,A 2 ). The 
minimal value of the t-channel momentum transfer is 
A?nin ~ — 4M^,?7 2 /(1 — rj 2 ). The unique feature of the 
process is that GPDs can be studied as functions of all 
three arguments independently. Note that if the timelike 
virtuality is away from the resonance region, the hadronic 
component of the photon can be neglected and conse- 
quently all subprocesses with vector meson production. 
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FIG. 3: The £-?7 coverage with E = 11 GeV electron beam. 
Only the regime Q 2 > M?i is shown. The kinematics is taken 
as follows: A 2 = -0.3 GeV 2 , Q 2 varies from 1 GeV 2 to 4 
GeV 2 and M 2 a from to Q 2 - lGeV 2 . The perimeters of 
areas covered for different lepton energy losses y are shown as 
described in the legend. 

As Fig. |^ demonstrates, the region of the GPD surface 
accessible in this reaction for a fixed target experiment 
setup is quite extensive. For a given £, cf. Eq. (||), with 
xb = Q 2 1 (2M^Ey), the variable r\ varies in limits 

min^l < \r]\ < min{max|77|, |77 cut |} , (11) 

where rj is evaluated according to Eq. (||) and the bound 
|?7 cu t| = \f— A 2 /(AM^ — A 2 ) comes from the condition 
| A 2 1 > |A 2 lin |. The minimal value of £ is governed by the 
requirement of applicability of perturbative QCD treat- 
ment, Eq. (0), chosen here as Q 2 > M 2 g + lGeV 2 . The 
restriction |?7| > £ is a simple consequence of the timelike 
nature of the final state photon. 



The ratio of (|10|) to the DVCS signal 

1 / *- i/_\rr a c _ 

~3tT 
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" s , dAa 
dMl 



im'f 
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In 



Q 2 



(12) 



is of order 0.01 for Q 2 ~ 2 GeV 2 . Difficulties in mea- 
suring such a small cross section will be overcome at 
high-luminocity machines, like JLab@12GeV with L = 
10 35 cm _2 s _1 or in collider experiments due to growth of 
the cross section with increasing £. 
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FIG. 4: Azimuthal angle dependence of the scaled cross sec- 
tion of e~p — > e~pe~e + for spacelike (top) and timelike (bot- 
tom) q 2 at E — 11 GeV. The solid, dash-dotted and dashed, 
dotted curves represent predictions for FPD and DD (with 
6=1, _B sca — 9 GeV -2 ) models of Ref. 0, respectively. 



Top: Mf- t 



for solid, dashed and M, 



0.7 GeV 

^2 



for 

dashed, dash-dotted curves, respectively. Bottom: Qf" = 
for solid, dashed and Q 2 = 1 GeV 2 for dash-dotted, dotted 
curves, respectively. 

The cross section (^) being directly proportional to 
GPDs is extremely sensitive to their skewedness de- 
pendence. Such kinds of measurements can easily es- 
tablish the credibility of our current understanding of 
mesonic-like components of GPDs and thus constrain the 
region which is lacking in the establishment the total 
spin sum rule |l3|. In Fig. ||, we present estimates for 
JLab@12GeV kinematics using two different models of 
GPDs: A model without skewedness dependence, which 
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FIG. 5: The dependence of scaled cross section for FPD 
(solid) and DD (dashed) models at <j> = tv/2 for M^dAa (left) 
and v/ Q 2 dAa (right) for the same kinematical settings as in 
Fig. top and bottom panels, respectively. 
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FIG. 6: Same as in Fig. displayed for HERA kinematics 
with y = 0.45, A 2 = -0.1 GeV 2 , and Q 2 = 4 GeV 2 (left) and 



M£ = 4 GeV 2 (right). 



corresponds to the conventional parton density taken to 
be the same at all values of 77, the so-called forward par- 
ton distribution (FPD) model. And a model based on 
a more sophisticated construction |L8| which takes the 
aforementioned parton densities as input and leads to a 
nontrivial ^-dependence. It is called the double distri- 
bution (DD) model. For specific details we refer to Ref. 
fl2|| . The most prominent way to discriminate between 
these models is to study the or Q 2 dependence of 
the cross section ([[o|) as shown in Fig. ||. One will be 
easily able to distinguished between different behaviors 
by varying the timelike photon virtuality over a short in- 
terval below the p-meson threshold. As Fig. |] shows, the 
higher energy of HERA results into even higher sensitiv- 
ity to the skewedness dependence of GPDs. We like to 
mention that at HERA the skewedness dependence can 
already be studied in measurements of the unpolarizcd 
cross section, since the imaginary part of the Compton 
amplitude dominates over the real part. 

To conclude, the exclusive electroproduction of lepton 
pairs provides a unique opportunity to determine exhaus- 
tive information on the nucleon's internal structure by 
accessing inter-parton correlations. The underlying GPD 
can be mapped uniquely as a function of all its variables 



which encode dynamics in longitudinal and transverse 
spaces. Our analysis demonstrates a high sensitivity of 
the lepton helicity asymmetry to the skewedness depen- 
dence of GPDs. Similar conclusions apply to nucleon spin 
asymmetries which again extract the interference term of 
Bether-Heitler and Compton amplitude and separate yet 
another combination of functions in question. The angu- 
lar distribution of the hnal state lepton pairs has a very 
rich structure and its measurement will lead to an in- 
dispensable complementary information on GPDs. One 
might expect that such studies at collider energies will pin 
down the skewedness dependence of gluon GPDs in the 
small-£ region and thus reduce theoretical uncertainties 
in diffractive leptoproduction of vector mesons. Finally, 
let us point out that the difference between the space- 
and timelike region is perturbatively computable and so 
the onset of the light-cone dominance can be elucidated 
as well. 
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